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ABSTRACT 

We estimate the masses of elliptical galaxies out to five effective radii using planetary nebulae and globular 
clusters as tracers. A sample of 15 elliptical galaxies with a broad variation in mass is compiled from the 
literature. A distribution function-maximum likelihood analysis is used to estimate the overall potential slope, 
normalisation and velocity anisotropy of the tracers. We assume power-law profiles for the potential and 
tracer density and a constant velocity anisotropy. The derived potential power-law indices lie in between the 
isothermal and Keplerian regime and vary with mass: there is tentative evidence that the less massive galaxies 
have steeper potential profiles than the more massive galaxies. We use stellar mass-to-light ratios appropriate 
for either a Chabrier/KTG (Kroupa, Tout & Gilmore) or Salpeter initial mass function to disentangle the stellar 
and dark matter components. The fraction of dark matter within five effective radii increases with mass, in 
agreement with several other studies. We employ simple models to show that a combination of star formation 
efficiency and baryon extent are able to account for this trend. These models are in good agreement with both 
our measurements out to five effective radii and recent SLACS measurements within one effective radii when 
a universal Chabrier/KTG initial mass function is adopted. 



1. INTRODUCTION 

There is strong evidence for the presence of dark matter 
in spiral galaxies where the rotation curves of their extended 
cold gas discs remain flat out to large radii. However, ellip- 
tical galaxies are generally free of cold gas so they lack an 
equivalent to the HI rotation curves of spiral galaxies. Fur- 
thermore, the stellar component of ellipticals is dominated 
by random motions so their kinematics are more difficult to 
model and the results are bedevilled by the mass-anisotropy 
degeneracy. 

The confirmation of dark matter haloes surrounding 
elliptical galaxies has largely been confined to the brightest 
galaxie s using either the X-ray em i ssion of their hot gas 
(e.g lLoewenstein & White! 119991: |Q' Sullivan & Ponmari 
2004b iHumphrev et al.1 120061: Uohnson et alj|2009t iDas et al. 
201C) o r strong lensing techniques (e. g. | Treu & Koopmansl 

| Gavaz 



2004 



iRusin & Kochanekl 12005b iGavazzi et alJ 12007b 



Koopmans et al.1 12009b lAuger et all 120 1 OallFaure et al J 
201 lllLeier et al.ll2.01 ll) . Stellar dynamical studies from inte- 



grated light sp ectra can also be use d to estimate dynamical 
masses (e.g. iGerhard et al.1 l200lt ICappellari et al.1 12006b 
iThomas et al.ll2007b iTortora et al.ll2009l) . but such studies are 
generally limited to within a couple of effective radii, R e g. 

To study the outer reaches of elliptical galaxies (i.e. beyond 
~ 2 effective radii) requires distant tracers such as planetary 
nebulae (PNe) or globular clusters (GCs). In particular, PNe 
can be used to trace intermediate mass galaxies, whereas GCs 
and other mass probes are generally biased towards more mas- 
sive syst ems. To make use of t he larger mass range probed by 
the PNe. lDouglas et af] (120021) developed a specialised instru- 
ment - the PNe Spectrograph (PN.S) - to study the kinemat- 
ics of these tracers in elliptical galaxies. The early results 
of this pr oject suggested a dearth of dark matter in ordinary 
elliptic als dRomanowskv et ail 12003b . However, iDekel et all 
(120051) showed that a declining velocity dispersion profile is 
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consistent with a massive dark halo if the tracer anisotropy is 
radially biased. More recent work utilising PNe to trace the 
kinematics of intermediate mass ellipticals find that although 
there is some evidence for the presence of dark matter haloes, 
the fraction of dark matter (within 5 Re g) is somewhat lower 
than their higher mass analogues (e.g. [Douglas et al.l 120071 
i Napolitano et al. 1 12009b INapolitano et al J 120111) . In addition. 
iNapolitano et al.l (120091) find that the dark matter halo con- 
centrations (i.e C200 or c v j r ) of these intermediate mass ellipti- 
cals are lower than ACDM predictions. However, estimating 
the dark matter halo mass and concentration at the virial ra- 
dius requires a large extrapolation from t he radial range of 
the cu rrent data (5R ff ~ 0.1r v ; r ). In fact. iMamon & Lokasl 
(120051) caution that extrapolation of dynamical studies within 
~ 5i? ff to the virial radius are fraught with large uncertain- 
ties. A comparison between intermediate and high mass ellip- 
ticals is difficult as it is rare for both mass scales to be probed 
using the same method. 

Whilst the overall mass of a galaxy can be derived using 
dynamical modelling of tracer kinematics, one must disen- 
tangle the stellar component in order to study the properties 
of the dark matter halo. However, incomplete knowledge of 
the initial mass function (IMF) inhibits this decomposition. 
While studies of our own Galaxy favou r an IMF with a flat- 
tened slo pe be low O.5M (IScalolll986b iKroupa et all [1991 
iChabrierlEOOl . near infrared spectroscopic studies of mas- 
sive ellipticals find that the low mass slope o f the IMF may 
become steeper (Ivan Dokkum & Conrovll20Tol) . Furthermore, 
seve ral recent studies have suggested a m ass-dependent IMF 
(e.g. lAuger etllbOlObbfTreu et al.ll20Tob . 

Clearly, our knowledge of the dark matter haloes surround- 
ing elliptical galaxies is far from complete. In this work, 
we study the outer regions of elliptical galaxies over a range 
of masses. These regions are relatively unexplored, espe- 
cially for the less massive systems. To this end, we com- 
pile a sample of galaxies with kinematic tracers beyond ~ 
2 Rpff . Previous authors have used either Jeans m odelling 
(e.g. INapolitano et all 12009b INapolitano et all 120 111) or orbit 
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Table 1 

Parent galaxy properties - (1) Name of galaxy. (2) Type of tracer (GC or PNe). (3) Distance i n Mpc. (4) and (5) Right asc ension and declination. (6) Position 
angle. (7) Effective radius . (8) B band luminosit y. The B band luminosity is derived from the de Vaucouleurs et al. 1 1991) apparent magnitudes and corrected 
for extinction according to Schlegel et al. 1 1998). (9) Heliocentric velocity. (10) Single power-law density index of tracers. (1 1) KS probability of approximate 
single power-law begin drawn from observed density model (12) Number of tracers beyond two effective radii. (13) Fraction of kinetic energy in rot ation. (14) 
References. The s ource of the effective radi i and tra cer density distribution are indi cated by the f and * symbo ls respectively: Capoellari et al. 2006 (Cap()6), 
ICoccato et alJ2009l(Coc09 ).ICote et al. 2001 (CotOl). ldeTor enzi et al. 20081 (De 08).IDirsch et alJ2005l (Dir05) .IDouglas et all2007KD ou07). IForbes et alJ20Qi 
(F06) , IHwang et a l. 2008 (H08). ILee et aU2008l (L08).IMei et al.1 2007 (Mei()7), Mendez et al. 2009 (Men()9), Napolitano et al. 2009 (N09), Napolita no et al. 
201 1 (Nl 1), Pens et al. 2004 (P04), Romanowskv et al. 2009 (R09). Saalia et al. 2000 (SagOO), Schuberth et al. 2006 (Sch06), Schuberth et al. 2010 (SchlO) 
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Teodorescu et al. 2005 (T05), Teodorescu et al. 20 1 1 {Til ), Woodley et al. 2010 (W10) 



librar y techniques, such as Schwar zschild (e.g. iThqmas et al.1 
12011b or NMAGIC modelling (e.g lde Lorenzi et al.ll2009l) . to 
study the kinematics of such tracers. Whilst effective, the lat- 
ter methods are arduous and generally applied on a galaxy by 
galaxy basis. In the coming years, where the number of kine- 
matic tracers surrounding elliptical galaxies is likely to dra- 
matically increase, it is important to develop methods to anal- 
yse a large sample of systems both quickly and effectively. 
Herein, we adopt a distribution function analysis. The advan- 
tage of such a scheme is that it allows us to study a number of 
systems with relative ease so the overall trends with mass can 
start to be addressed. 

The paper is arranged as follows. In §2, we introduce our 
sample of elliptical galaxies with distant kinematic tracers 
compiled from the literature. §3 describes the distribution 
functions and maximum likelihood analysis. We give our re- 
sults in §4 and develop some simple model predictions. Fi- 
nally, we draw our main conclusions in §5. 

2. ELLIPTICAL GALAXY SAMPLE 

Our aim is to probe the dark matter haloes of early type 
galaxies. To this end, we construct a sample of local galaxies 
with tracers reaching beyond two effective radii. Our sam- 
ple of 15 galaxies is compiled from the literature - the prop- 
erties of these systems are given in Table Q] along with the 
associated references. This sample covers a range of galaxy 
masses and environments (i.e. from field to cluster galaxies). 
The tracers are either GCs or PNe. In two cases (NGC 4649 



and NGC 5128), we use both GCs and PNe as tracers. As 
the different tracers may have different dynamical properties 
(e.g. anisotropy), we analyse each sample separately. How- 
ever, we do not attempt to model red and blue globular cluster 
populations separately. Whilst previous authors have found 
these populations may have different density profiles and or- 
bital properties (e.g iCote et al.l |200TI iHwang et all 12008b 
ISchuberth et al.ll2010l) . we choose to study the globular clus- 
ter population as a whole to maximise the number of tracers. 
Many of the PNe samples derive from the PNe Spectrograph 
projecfl This project specifically targets the outer regions of 
local galaxies and promises to increase substantially the num- 
ber of systems with dynamical tracers in the near future. 

We exclude any obvious outliers in the samples using a 3-a 
velocity clipping method (see e.g. iDouglas et al.ll2007l) . The 
line of sight velocities of tracers at similar projected radii are 
used to exclude any objects with obviously inflated veloci- 
ties (i.e. those with velocities exceeding 3a). In addition, we 
exclude outliers flagged in the literature. For example, we 
exclude those PNe fr om the NGC 3379 sample identified by 
IDouglas et ail (120071) as belonging to NGC 3384. 

For each sample of tracers, we approximate the density pro- 
file by a single power-law p oc r~ a . We use the density pro- 
files given in the literature. Where 2D profiles are given, we fit 
a single power-law to the projected distribution and increase 
the power-law index by 1 to convert to a 3D spatial profile. 

2 http ://w w w. strw.leidenuniv.nl/pns/PNS .public. web/PN . S .proj ect.html 
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The approximate single power-law is only fitted between the 
range of radii we are probing (e.g. between 2i? c ff and i? max )- 
In Table Q] we give the Kolmogorov-Smirnov (KS) probabil- 
ity that our approximate single power-law is drawn from the 
density profiles given in the literature. This statistic is omit- 
ted if the literature profiles are given as single power-laws. In 
most cases, a single power-law is a good approximation (with 
probability ~ 1). However, in NGC 4564 and NGC 4486 the 
probability is only ~ 0.3. Inspection of these profiles shows 
that the disagreement is driven by the outer 10% of tracers and 
hence the approximate single power-law is a good represen- 
tation of 90% of the sampl^B- In this exercise, we use the de- 
rived density profiles in the literature rather than the published 
positional data. Many of these density models are derived by 
taking into account completeness and selection biases (e. 
iDirsch et alj2005HDouglas et al. 120071: 1 Schuberth et al.l201 
and this information is not as readily accessible as the pub- 
lished positional and velocity data. However, we have also 
performed KS tests relative to published tracer positions and 
find similar results. 

Most tracers have a density power-law in the range a = 
2 — 4. In general, the GCs have more extended distribu- 
tions than the PNe. This may be related to the parent galaxy 
properties as the GCs tend to trace the more massive galax- 
ies. However, this may also be due to differences between the 
GC and PNe populations themsel ves. In general, PNe tend to 
follow the light dist ribution (e.g. ICoccato eTa l. 2009, but see 
iDouglas et al.ll2.007l), bu t the GCs can often be more extended 
(e.g. IDirsch etal .112005b 

In Fig. [Tj we show the size-mass relation of our sample (as- 
suming a Chabrier/KTG IMF) and the relation between tracer 
power-law index and stellar mass. As can be seen in these 
panels, more massive galaxies tend to have more extended 
tracer populations and larger effective radii. This size-mass 
relation proves im port ant to our conclusions and is discussed 
further in Section 14.11 Note that the effective radii we have 
adopted (given in Table [TJ are generally consistent with de 
Vau couleurs it! 1 / 4 fits . This eases comparison with studies 
like iHyde & Bernardil (120091) who adopt de Vaucouleur sur- 
face brightness profiles. However, we caution that the effec- 
tive radii can change substantially if a general Sersic i? 1//rl 
profile is adopted. 

3. DISTRIBUTION FUNCTION 

We analyse the dynamical properties of the tracers using 
a distribution function (henceforth DF) method. DFs are a 
valuable tool for studying steady state-systems as they replace 
the impracticality of following individual orbits with a phase- 
space probability density function. We provide a brief de- 
sc ription of these DF s bel ow but direct t he int erested reader 
to IE vans et all d 19971) and iDeason et al.l d201 lb . where more 
detailed descriptions are given. 

For simplicity, we use power-law profiles for the tracer den- 
sity and potential, namely p oc r~ a and $ oc r~ 7 , where a 
and 7 are constants. Although our formulae hold good for 
7 e [—1,1], models with 7 < are less useful for modelling 
galaxies. The velocity distribution is given in terms of the 



binding energy E = $(r) 



and the total 



Note that an increase (or decrease) of the tracer density power-law index 
by 1 dex leads to an increase (or decrease) in the mass estimate of ~ 30%. On 
average, our mass estimates are known to 20% (see Section 4); hence, only if 
the power-law is changed by 1 dex (a gross overestimate of the uncertainty) 
can the masses increase or decrease by more than the statistical errors. 
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Figure 1. Top panel: Size-mass relationship for the elliptical galaxy sam - 
ple. The solid line shows the quadratic model of Hyde & Bernardi ( 2009) 
and the dashed lines indicate a scatter of 0.2dex per stellar mass bin. Here, 
a Chabrier/KTG IMF is adopted. Bottom panel: Tracer power-law density 
index vs. stellar mass. The red triangles and blue squares indicate GCs and 
PNe respectively. Higher mass galaxies tend to have tracers with shallower 
density profiles. 
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where 

f(E) = £/3(7-2)/ 7 +a/ 7 -3/2 ^) 

Here, (3 is the Binney anisotropy parameter (lBinneviri980l) . 
namely 

a 1 («i> + H) 



2(v*) 



(3) 



which is constant for the DFs of form eqn (Q~|i. 

Our analysis assumes spherical symmetry. While the ma- 
jority of our elliptical galaxies look spherical in projection 
(E0/E1) there are a non-negligible number which are more 
flattened (E6). However, we only consider tracers beyond 
2Res, well beyond the region from which the galaxy type 
was inferred. With no prior knowledge of the shape of the po- 
tential in these regions, we make the simplest assumption of 
spherical symmetry. Relaxing this assumption calls f or more 
sophisticated modelling (e.g. Ide Lorenzi et alJl2007h . How- 
ever, we note that it is not obvious whether this extra compli- 
cati on makes any appreciab le difference to the mass estimates 
(e.g lde Lorenzi et alj|2009l) . 

These distribution functions can easily be ad apted to probe 
the r otational properties of the tracers (see IDeason et alj 
12011b . However, in the spherical approximation, it is un- 
physical for the tracer populations to have substantial rota- 
tion. In the limit of mild or no rotation, the mass profile and 
anisotropy are unaffected by the odd part of the distribution 
function. To this end, we proceed under the assumption that 
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the tracers are dominated by random rather than systematic 
motion. In Table Q] we estimate the fractional kinetic en- 
ergy in rotation, (V m t/<r) 2 , from the kinematic information 
give n in the literature fo r each sample of tracers (e.g. Table 
7 in lCoccato et al. 2009). In most cases, this fraction is small 
(< 20%) so it is safe to ignore rotation in our analysis and 
assume spherical symmetry. However, in a few cases (NGC 
4564, NGC 4649 (GC) and NGC 5128 (PNe)), the rotation is 
quite significant and our assumptions may not be valid. This 
is most apparent for NGC 4564 where V ro t/a ~ 1. It is be- 
yond the scope of this paper to model the tracer populations 
with oblate or triaxial spatial distributions, but we note that 
excluding these systems from our sample has little difference 
on our main results. 

3.1. Line of sight velocity distribution 

In our sample of local galaxies, we do not possess full six- 
dimensional phase space information for each tracer. How- 
ever, we can easily marginalise over the unknown compo- 
nents using the DF. We marginalise over the tangential ve- 
locity components (vi and Vb) and line of sight distance (D) 
to obtain the line of sight velocity distribution (LOSVD): 

F(l, b, v ios ) = J J J F(l, 6, D,vi, v b ,v ios )dvi dv b dD 

(4) 

Here and b are the Galactocentric longitude and latitude, 
whilst v\ os is the line of sight velocity. We assume that all the 
tracers are bound. Thus, the marginalisation over the tangen- 
tial velocity components and line of sight distance requires 

"tot < v CBC . 

This method can be generalized to account for observa- 
tional errors in the line of sight velocities. We assume the 
errors are Gaussian and marginalise over the line of sight ve- 
locity 

1 f (v-vios,o) 2 

F(l,b,v los )= / F(l,b,vi oa )e 2CT '™>° dv los 

y/ 27T(Tlos,0 J 

(5) 

Here, vi os ,o is the predicted line of sight velocity and ui os is 
the associated error. There are four parameters in our anal- 
ysis: the potential normalisation ($o), the potential power- 
law slope (7), the velocity anisotropy parameter and the 
tracer density power-law slope (a). We set the tracer density 
from the power-law approximations made in the previous sec- 
tion. There remains three parameters which we find using a 
maximum likelihood method. The likelihood function is con- 
structed from the LOSVD 

N 

\ogL(P, $0, 7) = Yl lo sHk, b i> ' Wl ° s - ' P> $ o, 7), (6) 

i=l 

Eqn © gives the three dimensional likelihood as a function 
of f3, 7 and $o- The total mass within a given radius can then 
be found from these parameters 

Here, $ n = $(lk pc). The analytic mass estimators of 
iWatkins et ail j2010l) also assume power-law models for the 
potential and tracer density. Our analysis extends beyond this 
formalism as, in addition to an estimate of the total mass, 
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Figure 2. The maximum likelihood contour levels for NGC 1 399 and NGC 
1344. The gray shaded regions show the la (68%) confidence region whilst 
the black lines encompass the 2<r (95%) confidence region. The top right 
panel of each figure shows the likelihood as a function of total mass within 
5i? eff . 



we also constrain the slope of the potential and the velocity 
anisotropy of the tracers. This is an important improvement 
as these parameters are often poorly constrained both in the- 
ory and observations. Furthermore, by constraining the slope 
of the potential (7), we can report results at a common radius. 

4. RESULTS 

We summarise the maximum likelihood parameters in Ta- 
ble |2 In Fig. I2 we show the likelihood contours for NGC 
1399 and NGC 1344 as examples. The grey shaded region 
shows the 68% confidence boundary and the solid line gives 
the 95% confidence boundary. Note that each panel shows a 
2D slice of the likelihood values. There is a strong degener- 
acy between the potential normalisation and power-law slope 
(note the 'banana' shape in the bottom right-hand panels). 
While we do not strongly constrain these individual param- 
eters, the mass profile is better defined as shown in the top 
right hand panels of Fig. [2] 

First, we com pare our results to th e analytic mass estima- 
tors described in W atkins et al.l(l2010h . These mass estimators 
compute the mass within the maximum 3D radius of the trac- 
ers, f max . For a spherical distribution of stars with density 
distribution p, the average 3D radius of a star with projected 
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[1o.j,jU.VJ 


UJ 


m ft ft 11 


rA A A Al 


/I A OA 

44o0 


o.U 


1 1 
-1.1 


r 1 o ft 11 

l-i.y,-u.ij 


r y( /: a ii 


4z.j 


[29.2,48.1] 


[zo.U, /Z.4J 


r\ a 
U.4 


TA TH71 

[U.2,U. /J 


TA 1 A 01 

[U.l,U.oJ 


/I /I O/l 

44V4 


/.4 


r\ 1 
U.l 


[-U.2,U.4J 


r A c A ci 


5.3 


r 1 a 'j m 
[1.9,3.9] 


r 1 c z h~\ 
[1.0,J. /J 


u. / 


in £ i ai 
[U.o,l.UJ 


TA 1 1 HI 
[U.J,1.UJ 


4 j 04 


1 A C 

1U.J 


i i 
-1.1 


r t i a oi 
[-2. 1,U.2J 


r A A A ci 


1.0 


[ft Q 1 71 


[U. o,Z. OJ 


u. / 


rA H 1 Al 


rA 1 1 ai 
[U.1,1.UJ 


40j0 


(3.0 


U.Z 


[U.U,U. J\ 




1 J.o 


[9.5,17.8] 


[8.3,26.3] 


U.4 


rn i a ci 


rn a a 71 

[U.U,U. / J 


4649 (GC) 


5.3 


-0.7 


[-1.5,0.2] 


[-3.4,0.6] 


14.0 


[7.2,15.8] 


[6.2,31.6] 


0.5 


[0.1,0.7] 


[0.1,0.9] 


4649 (PNe) 


4.1 


0.1 


[-0.2,0.4] 


[-0.6,0.5] 


18.0 


[6.9,21.7] 


[6.5,33.5] 


0.6 


[0.6,1.0] 


[0.0,1.0] 


4697 


5.2 


-0.5 


[-1.0,-0.1] 


[-1.5,0.1] 


4.0 


[2.5,4.4] 


[2.3,6.0] 


0.7 


[0.6,1.0] 


[0.2,1.0] 


5128 (GC) 


5.4 


0.2 


[-0.1,0.4] 


[-0.4,0.5] 


9.2 


[6.2,10.4] 


[5.3,13.4] 


0.4 


[0.1,0.5] 


[0.0,0.8] 


5128 (PNe) 


15.6 


0.5 


[0.4,0.6] 


[0.3,0.7] 


10.7 


[6.1,12.6] 


[5.0,16.8] 


0.7 


[0.6,0.9] 


[0.4,1.0] 


5846 


3.9 


0.2 


[0.0,0.5] 


[-0.7,0.7] 


23.8 


[8.2,28.1] 


[6.8,53.2] 


0.6 


[0.6,1.0] 


[0.1,1.0] 



Table 2 

Likelihood parameters. We give the galaxy name and maximum projected radius (scaled by effective radius) in the first two columns. The remaining columns 
give the weighted mean of the velocity anisotropy, potential normalisation and potential power-law slope with their 68% and 95% confidence intervals. 



radius, R is given by: 



13.0 



ST^rpjR, z) dz 

I-oo P( R > Z ) dz 



(8) 



For a power-law distribution of tracers, this equation is ana- 
lytic and can be expressed as 



('•) 



r[a/2]r[a/2 - 1] 

(r>/2-i/2]) 2 



R 



(9) 



Here, a is the power-law index of t he tracer dens i ty pro - 
file. Using equations 26 and 27 from [Watkins et al.1 (120101) . 
we estimate the mass using the maximum likelihood 7 and 
j3 valuefl We also evaluate the mass assuming an isother- 
mal potential (7 = 0j3 or Keplerian potential (7 = 1) with 
isotropic orbits (/? = 0). These values are compared to our 
maximum likelihood mass within (r max ) (see eqn|7| in Fig. 

El 

We see ver y good agreem e nt bet ween our estimated masses 
and those of IWatkins et al.l (120101) when the potential slope 
and anisotropy are given. The agreement is within ~ 10% for 
all the systems. However, when an isothermal potential (blue) 
or Keplerian potential (red) with isotropic tracers is assumed, 
there can be quite large discrepancies. The mass is systemati- 
cally overestimated when an isothermal potential is assumed. 
In this case, the discrepancies can be as large as 100%. Con- 
versely, if we had adopted a Keplerian potential the mass 
would be systematically underestimated. The disagreement 



4 Note that Watkins et al. 12010) label the power-law potential index and 
tracer density index as a and 7 respectively. This is the opposite to the nota- 
tion adopted in this work. 

5 In the case 7 = 0, the potential is logarithmic (e.g. see equation 9 of 
Watkins et al. 2010). This is the limit of the power-law models as 7 — ¥ 0. 
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Figure 3. Comparison with the mass estimator of Watkins et al. 1 2010). We 
show the agreement between our masses within r m ax (M(r < r max )) and 
the estimated masses from the Watkins et al. 12010) formalism (M cs t). The 
black squares show the Watkins et al. (2010) mass estimates given our maxi- 
mum likelihood /3 and 7. The red triangles and blue diamonds show the mass 
estimators assuming velocity isotropy (J3 = 0) with a Keplerian potential 
(7 = 1) and an isothermal potential (7 = 0) respectively. The solid line 
shows a one-to-one relation, and the dashed lines show 10% discrepancies. 

is worsened for systems where the velocity anisotropy devi- 
ates from isotropy and/or the potential is not close to isother- 
mal or Keplerian. These findings illustrate the effectiveness 
of our technique to estimate masses, as we make no strong 
assumptions about the velocity anisotropy or potential power- 
law slope. 

In Fig. UJ we show the velocity anisotropy distributions of 
the globular cluster tracers (red lines) and PNe tracers (blue 
lines). Each measurement is described by a Gaussian centred 
on the estimated value with a dispersion given by the pre- 
dicted (la) error. These Gaussian kernels are then summed to 
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Figure 4. The velocity anisotropy distributions for the globular clusters (red 
dashed lines) and PNe (blue solid lines). The thick lines show the sum of the 
Gaussian kernels for each measurement. Individual contributions are shown 
with the thinner lines. In general, the PNe orbits are isotropic/mildly radially 
biased whilst the globular cluster orbits are isotropic/mildly tangential. 




M(r<5R e „)/M sol 

Figure 5. The slope of the overall power-law potential (<E> oc r -7 ) vs. the 
total mass enclosed within 5 effective radii. The dark and light gray bands 
show the 68% and 95% confidence regions respectively. The solid black, 
dot-dashed blue and dotted red lines give the model predictions (within 10 — 
lOOkpc) for bulge+halo, halo and bulge models respectively. 

produce an overall distribution which is shown by the thick 
lines. Individual contributions are shown with the thinner 
lines. The PNe are generally more radially biased than the 
globular clusters. This has been noted by previous authors 
who find that GCs generally have mildly tangential/isotropic 
orbits whilst PNe tend t o have mildly radial/isotropic or- 
bits (e.g. iHwang et alJ 120081: iRomanowskv et al] l2009t 



fidence regions respectively. The solid and dashed lines give 
model predictions for the power-law slope with their associ- 
ated sc atter. These models are described in more detail in 
Section 14.1.11 The black lines show the slope-mass relation 
for bulge+halo models where a Hernquist profile is assumed 
for the stellar component and a Navarro-Frenk-White profile 
is assumed for the dark matter component. The red and blue 
lines show the relation with only a stellar or dark matter com- 
ponent respectively. The slopes for these models are com- 
puted between 10 — 100 kpc. 

At present, our 7 values are too poorly constrained to dis- 
tinguish between different model profiles. In the 95% confi- 
dence interval the slopes lie in between an isothermal (7 = 0) 
and Keplerian (7 = 1) regime. However, there is tentative ev- 
idence for a trend with galaxy mass - namely, the less massive 
galaxies have steeper potential profiles than the more massive 
galaxies. This trend is clear in the models, but is less so in 
the data. We tested the (anti)-correlation of 7 with mass in the 
data using a Spearman rank test. Given the probability distri- 
bution of 7 values (derived from the marginalised likelihood 
distributions), we draw a 7 value at random for each galaxy. 
The correlation with mass is then tested with a Spearman rank 
statistic. The exercise is repeated for 10 5 trials. We find that 
only 20% of the trials have a 95% significant trend. This exer- 
cise suggests that it is premature to claim an (anti)-correlation 
between pot ential slope and mass from the current constraints. 
We note that lBarnabe et al.l (12.0 1 lb also find a tentative indica- 
ti on that less massive g alaxi es have steepr po t ential s lopes. 

lAuger et all d2010ab and iKoopmans et al.l (120091) find al- 
most isothermal profiles for the SLACS elliptical galaxy sam- 
ple within li? c ff- The SLACS sample is biased towards the 
more massive galaxies where isothermal models are also pre- 
dicted by the bulge+halo models. 

4. 1 . Dark matter fractions 

The maximum likelihood normalisation (i.e. total mass) 
and power-law slope define the overall potential. To sepa- 
rate the dark matter and baryonic co mponents, we as sume 
the stars follow a Hernquist profile dHernquistl 1 1 9901) . To 
convert lumino sity into stellar mass, we assum e either a 
Chabrier/KTG dChabried2003t iKroupa et d][l993b IMF with 
{Ml Lb)* = 3.5 - 5.5 or a Salpeter dSalpeterl 19551) IMF with 
{M/L bY =6-10 (e.g. iGerhard et alJl2001b iTortoraeFaT] 
I2Q09PL We assume the stellar mass-to-light ratios cover the 
given range with a flat prior, so all values in the range have 
equal weight. This range is taken into account in the esti- 
mated error of the stellar mass (see Table [3}. The fraction of 
dark matter within a certain radius is given by 



/dm(< r) = 1 - M*(< r)/A/ tot (< r). 



(10) 



I Woodley et al.l 120 1 Ot iNapolitano et all 120 1 ll) . Note that the 
most tangentially biased PNe system belongs to NGC 4564 
which, as noted in the previous section, has evidence for sub- 
stantial rotation. 

The slope of the overall power-law potential (<& oc r~ 7 ) is 
shown in Fig. [5] as a function of total mass within five effec- 
tive radifl The black points give the weighted mean 7 valuefl 
and the dark and light gray bands show the 68% and 95% con- 

6 The total mass is computed within a 3D distance of 5R c g . The scaling 
relation between 3D deprojected half-light radius (Tim) and 2D projected 

half-light radius is ~ 1.3 (e.g. Ciotti 1991). Hence, a spherical distance of limited to within one effective radius. Here, we concentrate 

5-f? c ff corresponds to m 3.8r]y 2 - 



Here, we have have assumed that the gas mass is negligi- 
ble relative to the stellar mass. As we are considering ra- 
dial scales of ~ O-lf yir this is a good approximation (e.g. 
IP' Sullivan et"aT]l2007l) . 

In recent years, the fraction of dark matter as a function of 
galaxy mass has been stud i ed extensively in the l iterature (e.g. 
Padmanab han et all 120041: ICardone et al.ll2009t iTortora et al] 
2009t iCardone & TortoralboiOt lAuger et al J |201 Oat iGrillol 
20101: [Barnabe et al. 201 1J). Most of these studies have been 



7 The 'weighted' mean is given by, x 
given by the likelihood values. 



where the weights are 



8 Note that stellar masses are larger by a factor of r 
IMF is adopted rather than a Chabrier/KTG IMF 



1.8 when a Salpeter 
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Name M(< 5R cS ) 


M/Lb 


IMr 


71 r* 

M 


/DM 


(JNuC) 


r i nil i\/f 1 

Liu m q \ 


(< 5R cS ) 




r 1 a 1 1 a j 1 

[10 li M Q J 


(< 5i? c ff) 


0821 


2.3 ± 0.6 


O ± Z 


L 


1 O _l_ A 1 

i.z ± o.z 


a cn I A 1 1 

0.59 ± 0.1 1 








s 


2.1 zb 0.4 


A T7 1 A 1A 

0.27 ± 0.20 


1344 


2.6 ± 0.5 


1 A 1 n 

14 ± Z 


/-I 


0.8 ± U.Z 


A H A _l_ A AC 

U. /4 ± O.Oj 








c 

s 


1 C I A O 

1.5 ± 0.3 


A C A 1 A 1 A 

0.54 ± 0.10 


1399 


12.5 ± 1.8 


If 1 c 

35 ± 5 


L 


1 I AT 

1.6 ± 0.3 


A AA 1 A HI 

0.90 ± 0.02 








S 


2.8 ± 0.5 


A 1 A AO 

0.82 ± 0.03 


1407 


9.4 ± 1.3 


iy ± z 


c 


T 1 _l_ A A 

z. 1 ± 0.4 


A Ol _l_ A AO 

U.oZ ± 0.03 








s 


3.8 ± 0.7 


0.67 ± 0.05 


3377 


0.7 ± 0.2 


1Z ziz Z 


c 


no i a a 
0.3 ± 0.0 


A 1f\ _l_ A A*7 
0. /0 ± 0.0/ 








s 


0.5 ±0.1 


A A £L 1 A 1 O 

0.46 ±0.13 


3379 


1.3 ± 0.2 


9 ± 1 


c 


A £ I A 1 

0.6 ±0.1 


A C 1 I A A*7 

0.61 ± 0.07 








s 


1.1 ± 0.2 


A O 1 I A 1 ^ 

0.31 ± 0.12 


4374 


15.9 ± 1.9 


T 1 1 O 

31 ± 3 




Z.Z ± 0.4 


A OA _l_ A A 1 
0.89 ± 0.01 








s 


4.0 ± 0.7 


A OA 1 A A^ 

0.80 ± 0.02 


A A O/C 

4486 


31.7 ± 3.2 


4j ± 4 




J.J ± U.6 


A m _l_ A A1 

U.VZ ± U.U1 










C O 1 1 A 

5.8 ± 1.0 


A DC I A AO 

0.85 ± 0.02 


4494 


1.2 ± 0.2 


g i 1 
J x 1 




1 A 1 A ^ 

l.U ± U.z 


A n inn 
U.3Z zb U.1Z 








s 


1.8 ± 0.3 


-0.21 ± 0.22 


4564 


0.4 ±0.1 


/ ± Z 




A O _l_ A A 

U.z ± U.U 


A C A _1_ A 1 *7 
0.54 ± 0.1 / 








s 


0.4 ±0.1 


A 1 O I A 1 A 

0.18 ± 0.30 


4636 


10.8 ± 1.9 


A A _l_ a 

40 ± 6 


/-I 


1 O _l_ A O 

l.z ± O.z 


A A 1 I A AO 

0.91 ± o.oz 








c 

s 


1 I A /I 

2.1 ± 0.4 


A O A I A AT 

0.84 ± 0.03 


4649 


8.8 ± 1.3 


14 ± z 




G _l_ AC 

z.o ± U.J 


A *7/1 -4- A A/1 

U. /4 ± U.U4 








s 


4.9 ±0.9 


0.54 ± 0.07 


4697 


1.4 ±0.2 


7 ± 1 


c 


0.8 ± 0.2 


0.51 ±0.08 








s 


1.5 ±0.3 


0.12 ±0.13 


5128 


4.7 ±0.5 


17 ± 1 


c 


1.2 ±0.2 


0.79 ± 0.02 








s 


2.2 ±0.4 


0.63 ± 0.04 


5846 


11.3 ±2.7 


28 ±6 


c 


1.8 ±0.3 


0.87 ± 0.03 








s 


3.2 ±0.6 


0.77 ± 0.05 



Table 3 

Dark matter fraction parameters. We give the galaxy name, total 
mass-to-light ratio within 5R c g , adopted IMF (Chabrier/KTG (C) or 
Salpeter (S)), stellar mass and fraction of dark matter within 5-R c ff . 

on the fraction of dark matter within five effective radii - a 
region relatively unexplored in local elliptical galaxies. 

4.1.1. Models 

Before proceeding, we construct some simple models to 
show the expected relation between dark matter fraction 
(within a certain radius) and galaxy mass. We apply the fol- 
lowing steps: 

• Abundance matching is used to relate the halo mass 
(Af2oo) t o the stellar mass. W e use the prescription 
given by iBehroozi et al.l (120101) . The stellar mass-halo 
mass relation is given by their Equation (21) with pa- 
rameters listed in Table (2). These authors adopt a 
ACDM cosmology using the WMAP5 parameters. The 
scatter in stellar mass for a given halo mass is ~ 0.15 
dex. This abundance matching is applicable for a 
Chabrier/KTG IMF. To convert to a Salpeter IMF, the 
stellar masses are increased by a factor of 0.25 dex. 

• An NFW profile is adopted to model the density pro- 
file of the halo. This profile is fully described by two 



parameters: the halo mass and concentration. We re- 
late the con centration to the halo mass using the rela- 
tion given in lMaccio et al.l ( 12008b for WMAP5. At this 
stage, we have a fully defined dark matter profile for a 
particular stellar mass. The scatter in this relation for a 
given halo mass is ~ U.l dex. 

• Using the quadratic size-mass relation given in 
iHvde & Bernardil (120091) based on SDSS data, we cal- 
culate the effective radius of a galaxy from the stellar 
mass. The scatter for a given stellar mass is ~ 0.2 dex. 
The form of this size-mass relation is shown in the top 
panel of Fig, Q] against the observed values in our sam- 
ple. 

• Finally, we adopt a Hernquist profile for the stellar pro- 
file which is defined from the total stellar mass and ef- 
fective radius. 

We now have a dark matter and stellar profile defined for 
any particular stellar mass. In addition, we take into account 
the scatter introduced in each step of the analysis. These are 
propagated forward using Monte Carlo techniques. We em- 
phasize that there are no free parameters in these models that 
need to be fitted to the data. 

The models predict an increase in dark matter fraction 
within 5 i? c ff with increasing total/stellar mass (see black line 
in Fig. |5J. This general trend is independent of the IMF used. 
The fraction of dark matter within a scaled number of effec- 
tive radii is driven by two factors: the star formation efficiency 
and the concentration of the se stars within the dark matter 
halo (cf. IZaritskv et ai]|2008l) . 

The star formation efficiency is implemented in the mod- 
els from the stellar mass-halo mass relation. The abundance 
matching, by construction, reproduces the observable stellar 
mass function. It is well known that there is a universal U- 
shape d trend of star formation effici ency (e.g. | Bensonetal] 
2000t iMarinoni & Hudsonl I2002t iNapolitano et all 12005b 
van den Bosch et all 120071: IConrov & Wechsleri 120091) . The 
peak efficiency occurs at M* ~ 10 11 where increasingly mas- 
sive galaxies have a lower efficiency due to the large cooling 
time of their hot gas (e.g. IWhite & Reeslll978l) and the least 
massive systems are unable to retain their primordial gas con- 
tent for long enough to form stars. Hence, the the lowest mass 
and highest mass galaxies are the most dark matter dominated. 

The second factor driving the trend of dark matter fractions 
is the size-mass relationship of early type galaxies. With in- 
creasing total (or stellar) mass the effective radius increases 
(see Fig. [T]). In addition, the size increases more steeply 
with increasing mass. This 'baryon un-packing' is an im- 
portant factor when we are considering the dark matter con- 
tent within a scaled number of effective radii. Regardless 
of the dark matter halo properties (or star formation effi- 
ciency), a higher concentration of baryons (or smaller ef- 
fective radius) leads to smaller dark matter fractions. This 
simple scali ng behaviour has been no t ed by previous au- 
thors (e.g. iPadmanabhan et all 120041 iTortora et alj l2009t 
INapolitano et al.l l2010t iDutton et all 1201 lb . Note that disc 
galaxies have a much shallower size-mass relation. Thus, the 
dark matter fractions are approximately constant with mass 
for these galaxies (e.g. IDutton et al.112.01 lh . 

In Fig.|6l we illustrate the dependence of our model on the 
adopted stellar mass-halo mass relation, size-mass relation 
and cosmology. The solid black lines indicate our adopted 




Figure 6. The dark matter fraction within one (left panel) and five (middle panel) effective radii as a function of stellar mass. The solid black line is our adopted 
model. The dotted red line shows the relation when a linear stellar mass-halo mass relation is used M200 = AM*. The dashed blue line shows a model 
with a size-mass relationship applicable to disc galaxies. The dot-dashed green line shows the model prediction when a WMAP1 cosmology is adopted. In the 
right-hand panel we show the dark matter fraction as a function of radius (scaled by R c g) for the different models. The triangle and square symbols show the 
profiles for lower (M* ~ 2 X 10 10 Afo) and higher (M* ~ 5 X 10 JWq) mass galaxies respectively. 
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Figure 7. The dark matter fraction ratio within 5 effective radii vs. total 
mass within 5 effective radii (left) and stellar mass (right). A Chabrier/KTG 
((M/L)* ~ 4.5) and Salpeter ((M/L)* ~ 8) IMF is assumed in the top 
and bottom panels respectively. The solid line indicates the predicted relation 
from simulations. The predicted scatter is given by the dashed lines. The inset 
panels show the residuals of the data vs. models. The red square highlights 
NGC 4494 which is inconsistent with stellar mass-to-light ratios larger than 
(M/L)* ~ 5. 

model. We consider three variations to our adopted model 
(shown by the coloured dotted, dashed and dot-dashed lines): 



• Stellar mass-halo mass relation: We adopt a linear re- 
lation of the form A/200 = AM*. The normalisation 
constant is set using the mean halo mass corresponding 
to a stellar mass of M* ~ 10 11 from our original abun- 
dance matching relation. Models adopting this relation 
are shown with the red lines. 

• Size-mass relation: We consider a size-mass relation 
applicable to disc galaxies. This is derived using the 



iPizagno et al.l (120071) sample of disc galaxies. The stel- 
lar masses are de rived from the co lours according to 
the prescription of iBell et all (120031) . A -0. 15dex offset 
was applied to provide stellar masses consistent with 
a Chabrier/KTG IMF. We then fit a linear relation be- 
tween effective radius and stellar mass to approximate a 
size-mass relationship. Models assuming this size-mass 
relation are shown by the blue lines. 

• Cosmology: We adopt a WMAP1 cosmology instead of 
the more recent WMAP5 parameters. In th is case the 
stellar mass-halo mass relation is taken from lGuo et al.l 
(120101) (who adopt a WMAP1 cosmology) and the 
WMAP1 mass con centration relation is used from 
iMaccio et al.l (120081) . This is shown by the green lines. 

The purpose of these variations is to emphasize the key in- 
gredients in our adopted model driving the apparent trend. 
Thus, these variations are for illustrative purposes only and 
we are not proposing that they are viable alternatives to our 
adopted model. For example, adopting a size-mass relation 
applicable to disc galaxies is a particularly poor assumption 
when applied to elliptical galaxies. 

The left and middle panels show the dark matter fraction 
within one and five effective radii as a function of stellar 
mass. We have not included the scatter, so this figure only 
illustrates the mean dependence of the models. Note that a 
Chabrier/KTG IMF is assumed but the same trends are seen 
with a Salpeter IMF. Adopting a different cosmology only 
slightly affects the dark matter fraction-stellar mass relation. 
As alluded to earlier, the relationship is driven by star forma- 
tion efficiency and baryon extent. By adopting a shallower 
size-mass relation (applicable to disc galaxies) the dark mat- 
ter fractions are constant over a large range of stellar mass. 
A similar effect is seen by allowing for a linear stellar mass- 
halo mass relation. The effects are the most severe within one 
effective radii as the baryons are more dominant at smaller 
radii. 

In the right-hand panel, we show the dark matter fraction as 
a function of radius (scaled by i? c g ). The profiles of the higher 
mass galaxies (M* ~ 5 x 10 11 M Q ) are more model depen- 
dent than those of lower mass galaxies (M* ~2x 1O 1O M0). 

4.1.2. Comparison with data 
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In Fig. [7] we show the dark matter fractions within five ef- 
fective radii as a function of total mass within this radius (left) 
and as a function of total stellar mass (right). The top panels 
show the relation for a Chabrier/KTG IMF and bottom pan- 
els are for a Salpeter IMF. The symbols with error bars are 
the data points derived in this work and the solid and dashed 
lines give the predicted relation from the models with the as- 
sociated (la) scatter. The residuals of the data vs. model (i.e. 
(data-model)/model) are shown in the inset panels. 

In general, the model is in good agreement with the data 
when a Chabrier/KTG IMF is adopted. Our dark matter 
fractions are in good agreement with previous studies which 
find /dm(< 5i? c ff) ~ 0.4 — 0.5 for ordinary ellipticals 
(e.g. NGC 4494 NG C 3379, NGC 4697 see Figure 12 in 
iNapolitano etaill20TTl) and /dm(< 5i? cff ) ~ 0.8 - 0.9 for 
group or cluster central ellipticals (e.g. NGC 1399, NGC 
1407, NGC 44 86, NGC 4648, NGC 5846 see Figure 8 in 
iDas et al.ll20"Iol) . 

Note that our total mass estimates (given within 5i? c g in Ta- 
ble [3}, and hence dark matter fractions, usually have smaller 
uncertainties than t hose found in more general dynami- 
cal m odelling (e.g. Ide Lorenzi et al.l 120071: Ide Lorenzi et al.l 
120091) . This is because such numerical modelling proce- 
dures allow for more general solutions (e.g. non-constant 
anisotropy, triaxiality) and thus cover a larger parameter 
space. Nonetheless, our simplistic approach is much more 
easily applied to a large sample of galaxies and hence this 
work is complementary to more general numerical methods. 

The red square highlighted is NG C 4494, which has a cu- 
riously low dark matter fraction. Napolit ano et al.l (12009b 
studied this system in detail and found an abnormally low 
dark matter concentration and also derive a similar dark mat- 
ter fraction within 5i? c g. The remaining low mass galaxies 
have dark matter fractions consistent with the model predic- 
tions (at least for a Chabrier/KTG IMF). While the case of 
NGC 4494 is in contradiction to the predictions from simula- 
tions, the general trend suggests that the dark matter proper- 
ties of lower mass galaxies are not very different from those 
of higher masses. However, a larger sample of tracers, prob- 
ing both low and high masses, is needed to investigate this 
further. 

Adopting a Salpeter IMF causes some of the low mass sys- 
tems to have lower than predicted dark matter fractions. On 
the other hand, a Chabrier/KTG IMF is able to describe the 
trend over a wide range of masses. A secondary factor that 
may influence the dark matter fraction-mass relation is varia- 
tions in stellar population properties, such as a non-universal 
IMF (e.g lAuger et al.ll2010UlTreu et al.l2010l) . We find no ev- 
idence for a variation in IMF with mass from Fig. [7] The sug- 
gested IMF variation, where the stellar-mass to light ratio in- 
creases with mass, would require the fraction of dark matter to 
increase more steeply when a universal IMF is adopted. How- 
ever, we are probing out to five effective radii, well beyond 
the regions where the baryonic material is dominant so are 
less sensitive to IMF variations than studies probing within 
one effective radii. We compare our models with such studies 
in the following section. 

4.1.3. Comparison with SLACS data 

In Fig. [8] we show the projected dark matter fractions as 
a function of galaxy mass and ste llar mass. The data points 
derive from the SLACS sample of lAuger efafl (I2010al) . This 
work focuses on relatively high mass galaxies within one ef- 
fective radius. The simple models that we employ here are 
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Figure 8. The projected dark matter fraction ratio within 0.5 effective radii 
vs. total mass within 0.5 effective radii (left) and stellar mass (right). The 
data points are for the SLACS sample of elliptical galaxies published in 
lAuger et alj I2010af) . The inset panels show the residuals of the data vs. 
model. 

able to explain the overall trend of increasing dark matter 
fraction with galaxy mass. Similar to the case at larger radii, 
adopting a Salpeter IMF leads to lower dark matter fractions 
than the model predictions. In fact, several systems are con- 
sistent with a negative dark matter fraction, which suggests 
that adopting a Salpeter IMF for such systems is unphysical. 
Once again a universal Chabrier/KTG IMF is able to explain 
the dark matter fraction-mass relation reasonably well. Vari- 
ations in IMF with mass can in principle be probed by devia- 
tions of the data from the model predictions. Qualitatively, we 
see no evidence for a difference in slope between the data and 
the model. However, a more sophisticated model is required 
to investigate these variations in detail. In summary, we find 
that the variation in dark matter fraction with galaxy mass is 
consistent w ith a simple, u niversal IMF halo model. 

Recently, iGrillol d2010l) finds a constant, if not decreas- 
ing, projected dark matter fraction with galaxy mass. This 
is inconsistent with both our models a nd our meas u rement s 
as well as with the measurements by lAuger et alj d2010al) . 
independently of the adopted IMF. IGrillol d2010l) suggests 
that these contrasting results may be due to other studies be- 
ing more model de pendent (e.g. Padmana bhan et ail l2004t 
iTortoraet alj 12009 who give three dimensional dark matter 
fractions) and/or the more massive early type galaxy bias of 
their sam ple. However, this is difficult to reconcile with the 
results of lAuger et alj ( I2010al) who also give projected dark 
matter fractions (i.e. with minimal modelling assumptions) 
and have a similar mass bias in their sample. 

5. CONCLUSIONS 

We studied the mass profiles of local elliptical galaxies us- 
ing planetary nebulae (PNe) and globular clusters (GCs) as 
distant tracers (with R > 2R c g). A sample of 15 galaxies 
was compiled from the literature. A distribution function- 
maximum likelihood method was used to study the dynamics 
of the tracers under the assumptions of spherical symmetry 
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and power-law models for the (overall) potential and tracer 
density. We summarise our conclusions as follows: 

(1) We compared our distribution function-m aximum likeli- 
hood m ethod to the analytic mass estimators of IWatkins et al.l 
( 120101) . There is good agreement when we adopt our maxi- 
mum likelihoo d velocity anisotropy and potential power-law 
indices into the W atkins et al.l ( 120101) formalism. However, as- 
suming isotropy and an isothermal/Keplerian potential leads 
to a systematic overestimate/underestimate of the mass (by up 
to ~ 100%). Our method makes no strong assumptions about 
the values of the velocity anisotrop>0 and potential power-law 
slope. This is an important improvement, as these parameters 
are difficult to constrain by other means. 

(2) The PNe are generally more central ly concentrated (fol- 
lowing closely the stellar luminosity, see lCoccato et al.ir2009l) 
and are on more radial orbits than the GCs. This may be re- 
lated to the parent galaxy properties as the PNe tend to trace 
the less massive systems, but could also reflect differences 
between the GC and PNe populations (e.g. GCs are generally 
older). 

(3) The slope of the overall (power-law) potential is not 
strongly constrained but lies in between the isothermal (7 = 
0) and Keplerian (7 = 1) regimes. In the models the less 
massive galaxies have steeper potential profiles than those of 
the more massive galaxies. There is tentative evidence that 
this trend is also present in the data. However, a Monte 
Carlo test using the Spearman rank statistic shows no evi- 
dence for a clear tre nd. Recently. iKoopmans et al.l (120091) and 
lAuger et"aT] d2010al) found roughly isothermal density profiles 
for their sample of (massive) elliptical galaxies within an ef- 
fective radius. This is in good agreement with bulge+halo 
models for the more massive elliptical galaxies. 

(4) We constructed simple halo models to predict the expected 
relation between dark matter fraction and halo mass. The 
combination of star formation efficiency (stellar mass-halo 
mass relation) and baryon concentration (size-mass relation) 
is able to describe the observed trend of increasing dark matter 
fraction within a scaled number of effective radii. We find no 
evidence for an additional factor, such as a non-universal IMF, 
affecting this relation. The dark matter fractions are consis- 
tent with the models when a universal Chabrier/KTG IMF is 
assumed for elliptical galaxies. In particular, a Salpeter IMF 
is inconsistent with so me of the lower mass gal axies. This is 
in good agreement with lCappellari et ail (120061) who reached 
the same conclusion using a completely different approach. 

Finally, let us remark that we have neglected the influence of 
baryonic processes on the distribution of dark matter. Our 
adopted NFW profile is applicable for a 'pristine' dark matter 
distribution, un-modified by the process of galaxy formation. 
For example, collapsing gas can exert a gravitational drag on 
the dark matter leading to a more concentrated dar k matter 
halo (i.e. iBlumenthal et1Illl986tlGnedin et al.ll2004l) . On the 
other hand, rapid supernova driven feedback processes can 
drive the da rk matter particles away fro m the centre of the 
galaxy (e.g. iPontzen & Governatoll201 lb leading to halo ex- 
pansion. Halo contraction or expansion can increase or de- 
crease the fraction of dark matter within a given radius. The 

9 However, we do assume that /3 is constant with radius 



importance of these processes is poorly understood and we 
make not attempt to model them in this work. However, we 
note that this is another secondary effect (in addition to IMF 
variations) that should be explored with more sophisticated 
modelling. 
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